arXiv:1502.07066v3 [cs.IT] 21 Oct 2015 


1 


To Harvest and Jam: A Paradigm of Self-Sustaining 
Friendly Jammers for Secure AF Relaying 

Hong Xing, Kai-Kit Wong, Zheng Chu, and Arumugam Nallanathan 


Abstract —This paper studies the use of multi-antenna harvest- 
and-jain (HJ) helpers in a multi-antenna amplify-and-forward 
(AF) relay wiretap channel assuming that the direct link between 
the source and destination is broken. Our objective is to maximize 
the secrecy rate at the destination subject to the transmit power 
constraints of the AF relay and the HJ helpers. In the case of 
perfect channel state information (CSI), the joint optimization 
of the artificial noise (AN) covariance matrix for cooperative 
jamming and the AF beamforming matrix is studied using semi- 
definite relaxation (SDR) which is tight, while suboptimal solu¬ 
tions are also devised with lower complexity. For the imperfect 
CSI case, we provide the equivalent reformulation of the worst- 
case robust optimization to maximize the minimum achievable 
secrecy rate. Inspired by the optimal solution to the case of 
perfect CSI, a suboptimal robust scheme is proposed striking 
a good tradeoff between complexity and performance. Finally, 
numerical results for various settings are provided to evaluate 
the proposed schemes. 

Index Terms —Harvest-and-jam, artificial noise, cooperative 
jamming, amplify-and-forward relay, wireless energy transfer, 
physical-layer security, semi-definite relaxation, robust optimiza¬ 
tion. 

I. Introduction 

The pressing demand for high data rate in wireless commu¬ 
nications networks coupled with the fact that mobile devices 
are physically small and power-limited by batteries, has driven 
the notion of energy harvesting (EH) to become a promising 
resolution for green communications [1,2]. Among the var¬ 
ied available resources for EH, radio-frequency (RE)-enabled 
wireless energy transfer (WET) has aroused an upsurge of 
interest for its long operation range, ubiquitous existence in 
the electromagnetic radiation, and effective energy multicas¬ 
ting, which motivates the paradigm of simultaneous wireless 
information and power transfer (SWIFT), e.g., [3-6]. 

A typical SWIFT system consists of one access point (AF) 
that has constant power supply and broadcasts wireless signals 
to a group of user terminals, amongst which some intend 
to decode information, referred to as information receivers 
(IRs), while others scavenge energy from the ambient radio 
signals, named energy receivers (ERs). This gives rise to a 
challenging physical (FHY)-layer security issue where the 
ERs may eavesdrop the information sent to the IRs due to 
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their close proximity to the AF To overcome this problem, 
in [7-9], several researchers presented various approaches to 
guarantee secret communication to the IRs and maximize the 
energy simultaneously transferred to the ERs or to satisfy 
the individual EH requirement for the ERs and maximize the 
secrecy rate for the IR, by advocating the dual use of the 
artificial noise (AN) or jamming. 

However, previous works all assumed that the ERs in the 
SWIFT systems attempt to intercept the information for the 
IR, which is overly protective. On the contrary, it is possible 
that some ERs are cooperative, especially when they are EH- 
enabled wirelessly. Eollowing the recent advances in wireless 
powered communications networks [10,11], this paper pro¬ 
poses a self-sustaining harvest-and-jam (HJ) relaying protocol, 
where in the first transmission phase a single-antenna trans¬ 
mitter transfers confidential information to a multiple-antenna 
amplify-and-forward (AE) relay and power to a group of multi¬ 
antenna EH-enabled idle helpers simultaneously, while in the 
second phase, the relay amplifies and forwards the information 
to the IR under the protection of the AN generated by the 
helpers using the energy harvested from their received signals 
in the first transmission phase. 

Fhysical (FHY)-layer security issues in the rapidly growing 
cooperative networks have attracted much attention. Cooper¬ 
ative approaches, such as, cooperative jamming, communi¬ 
cations have been widely examined [12-15]. The idea is to 
assist the transmitter in the secrecy transmission by generating 
an AN to interfere with the eavesdropper via either multiple 
antennas or external trusted helpers [16-19]. However, all of 
those utilizing ANs require additional supply of power and 
therefore incur extra system costs. Meanwhile, collaborative 
use of relays to form effective beams jamming the eavesdrop¬ 
per, i.e., secure collaborative relay beamforming, has been 
studied for relay-wiretap channels with single eavesdropper 
in [20], multiple eavesdroppers with AE relays and decode- 
and-forward (DE) relays in [21] and [22], respectively. All, 
however, assumed the availability of perfect channel state 
information (CSI). Though [23] proposed robust AE relay 
beamforming against the eavesdropper’s channel, the solutions 
were yet suboptimal. 

The assumption of perfect CSI of the eavesdroppers appears 
to be too ideal because the eavesdroppers, despite being 
legitimate users, wish to hide from the transmitter without 
being cooperative in the stage of channel estimation. Even if 
they are registered users and bound to help the transmitter 
in obtaining their CSIs to facilitate their own communication, 
the CSIs at the transmitter side will change due to mobility 
and Doppler effect, and may be outdated. Moreover, even 
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for the legitimate users, the estimated CSIs may also be 
subject to quantization errors due to the limited capacity of 
the feedback channel, although the inaccuracy is reasonably 
assumed less severe than that for the eavesdroppers. To tackle 
this issue, state-of-art schemes have been developed ([24] 
and the references therein), among which the worst-case 
secrecy rate is commonly employed to formulate the robust 
secrecy rate maximization problem [8,19,25-27]. The robust 
transmit covariance design for the secrecy rate maximization 
in a multiple-input-single-output (MISO) channel overheard 
by multi-antenna eavesdroppers was considered in [25,28] 
while the enhanced secrecy performance was achieved by 
introducing a friendly jammer in the same scenario in [26], in 
which a joint optimization of the robust transmit covariance 
and power allocation between the source and the helper was 
studied via geometric programming. More recently, [8] studied 
a joint robust design of the information beams, the AN and 
the energy signals for SWIFT networks with quality-of-service 
(QoS) constraints. 

The contribution of this paper is threefold. First, with perfect 
CSI, in addition to the joint optimal solutions, we propose 
two near-optimal schemes with much reduced complexity by 
exploiting the optimal structure of the relay weight matrix, 
and providing a semi-closed form solution for the relay weight 
matrix given fixed null-space jamming, respectively. Second, 
besides the imperfect eavesdropper’s channel, legitimate chan¬ 
nels such as those from the K HI helpers (the transmitter) to 
the legitimate receiver {K HJ helpers), and from the AF relay 
to the receiver are jointly modeled with imperfect estimation, 
and multiple semi-indefinite non-convex constraints have been 
judiciously replaced by linear matrix inequalities (LMIs) to 
fit the semi-definite programming (SDP). Third, a rank-one 
reconstruction algorithm exploiting the structure of the semi- 
definite relaxation (SDR)-based solutions has been proposed 
to provide promising performance at low computational cost. 

Of particular relevance to our work is [29] which jointly 
optimizes the AF matrices and AN covariances in a relay 
wiretap channel with multiple multi-antenna AF relays and 
multiple multi-antenna eavesdroppers via a worst-case robust 
formulation. While our network model is similar, the differ¬ 
ence of our work from [29] is twofold. On one hand, in this 
paper, the AN generated by the friendly jammers are subject to 
their respective channels from the transmitter during WET in 
the first transmission phase. On the other hand, the technique 
in [29, Proposition 7] cannot be applied to our problem since 
the AN beams and the forwarded information are transmitted 
via different channels in ours. As a consequence, to the best 
of authors’ knowledge, our proposed worst-case based robust 
optimization scheme that incorporates imperfect CSIs into all 
the HJ helpers, has not been addressed in the literature. 

It is worth noting that devising a wireless-powered friendly 
jammer to enhance PHY-Iayer security for a direct transmis¬ 
sion protocol was studied in [30], in which the “harvesting” 
blocks and “jamming” blocks were well exploited to compose 
four different types of harvesting-jamming cycles. Compared 
to [30], which focused on the dedicated scheduling of “har¬ 
vest” and “jam” operations and its long-term performance, 
ours are concerned with adaptive rate/power optimization with 


multiple HJ helpers to achieve higher worst-case secrecy rate. 
Moreover, instead of assuming perfect channels to/from the 
HJ helpers, our robust optimization algorithm takes imperfect 
legitimate channels into account to provide robustness. 

Note that in this paper, as in [23,29], we assume that the 
channel between the transmitter and the AF relay is perfectly 
known and there is no direct link between the transmitter and 
the receiver or the eavesdropper, a common assumption in the 
concerned AF relay wiretap channel [20,21]. 

Notations —Throughout, we use the upper case boldface 
letters for matrices and lower case boldface letters for vectors. 
The superscripts (-j^, (-j^ and {■)^ represent the transpose, 
conjugate and conjugate transpose, respectively. Also, tr(-) 
and E[-] stand for the trace of a matrix and the statistical ex¬ 
pectation for random variables, respectively. Likewise, vec(A) 
is defined as a column vector obtained by stacking the rows 
of A on top of one another. vec^“^^ is the inverse operation 
of vec. null(A) denotes the null space of A. 0 represents the 
Kronecker product of two matrices. In addition, the notation 
A A 0 indicates that A is a positive semi-definite matrix and J 
(0) denotes an identity (all-zero) matrix with appropriate size. 
Furthermore, || • || represents the Euclidean norm of a vector, 
while Pr{-) stands for the probability of an input random 
event. Einally, [a;]+ denotes max(0,x) and (•)* stands for an 
optimal solution. 

II. Network Model 

We consider a cooperative relay wiretap channel for SWIFT 
over a given frequency band as shown in Fig. 1(a). We assume 
that there is a transmitter, named Alice, sending confidential 
messages to the IR, Bob, in the presence of an eavesdropper 
[31], Eve, with the aid of a multi-antenna AF relay and K 
ERs willing to act as HJ helpers, 77helper = 

The transmitter, ERs, and the AF relay are deployed in a same 
cluster that is relatively far away from the destination and Eve, 
such that there is no direct link from the transmitter to the 
receiver or Eve, respectively. Moreover, the ERs are assumed 
to be located closer to the transmitter than the AF relay in 
order that they can harvest sufficient amount of energy for 
jamming. Alice, Bob and Eve are all assumed to be equipped 
with single antenna, while the AF relay and each of the K 
helpers are assumed to have the same Nt antennas. 

Using two equal slots for the HJ relaying protocol, as shown 
in Fig. 1(b), for the first phase, Alice sends a confidential 
message to the relay while simultaneously transferring energy 
to the K helpers; for the second phase, the relay amplifies and 
forwards the message to Bob while the K helpers perform 
cooperative jamming using their respective harvested energy 
from the first transmission phase, to compromise Eve. In this 
paper, we assume a quasi-static fading environment and for 
convenience denote hg G C^* ^ ^ as the complex channel from 
the transmitter to the relay and hk G k = 1,... ,K, 

as that from the transmitter to the fcth helper; ho as the 
transpose of the complex channel from the relay to Bob 
and hk G k = as that from Hfc to Bob; 

go G and G fc = 1,..., 7T, as those from 

the relay and to Eve, respectively. 
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-^ First transmission phase - Energy transfer 

-► Second transmission phase X. _ ^ Jamming 



(a) AF-relaying wiretap channel with jamming. 



Fig. 1. HJ-enabled cooperative relaying for secure SWIFT. 

In the first transmission phase, the baseband received signal 
at the AF relay can be expressed as 

= ho^/^s + nr, (1) 

where s is a circularly symmetric complex Gaussian (CSCG) 
random variable, denoted by s ~ CAf{0, 1) and Ur is the 
additive complex noise vector, denoted by rir ^ CAf{0, cr^I)- 
Also, Ps denotes the given transmit power at Alice. Further, 
the received signal at each helper H^. is expressed as 

yk = hk\/~PsS + n't., ( 2 ) 

where is the additive noise, denoted by nj. ^ CA/^(0, crf^I). 

On the other hand, for WET, the harvested energy of in 
each unit slot is given by 

Ek = r,E[\\hk^/Ksf] = yPs \\hkf , Vfc, (3) 

where 0 < p < 1 denotes the EH efficiency. 

In the second transmission phase, the linear operation at the 
AE relay can be represented by 

x' = Wy,, (4) 

where x' G is the retransmit signal at the AE relay 

and W G is the beamforming matrix. Note that the 

transmit power of the AE relay can be shown as 

tr (E [xx^] ) = tr (W (Psh^h^ + , (5) 

which is constrained by the maximum available power at the 

AE relay, i.e., P^, which is given by 

tr (W (Pshoh^ + all) < Pr. (6) 

In the meantime, each Hfe will help generate an AN rik G 

QNtxi J.Q interfere with Eve. Similar to [16], we assume 


that rifc’s are independent CSCG vectors denoted by ^ 
CA/^(0, Qk), Vfc, since the worst-case noise for Eve is known 
to be Gaussian. In addition, each Hfc has a transmit power con¬ 
straint due to its harvested energy in the previous transmission 
phase, i.e., tr(Qfc) < yPgWhkW^ (c.f. (3)), Vfc. 

The received signal at Bob can thus be expressed as 

T ^ T T 

Vb = Whos + hk nk fig Wrir -l- Jif), (7) 

where nb ^ CAf{0,alI) is the additive noise at Bob. Simi¬ 
larly, the received signal at Eve can be expressed as 

K 

ye = \/^goWhos + Y9k'>^k+goWnr + ne, (8) 
fc=i 

where rie CAf{0, all). According to (7) and (8), the signal- 
to-interference-plus-nose-ratio (SINR) at Bob and Eve can be, 
respectively, expressed as 

- _ PslhWhol^ _ 

Jo •~T Tc' ~-T 

alK ww^hl + Eti K QkK + 

and 

. =__. (10) 

alglWW^gl -k Ylk=i alQkdl + 

As such, the achievable secrecy rate at Bob is [16] 

ro = ^[log2(l+7ft)-log2(l-k7e)]’^- (11) 

III. Joint AN-AE Beamforming with Perfect CSI 

A. Problem Formulation for Perfect CSI 

We aim to maximize the secrecy rate at Bob subject to the 
transmit power constraints at the AE relay and each individual 
helper H^, k = 1,..., K. Thus, our problem is to solve 

(PI) : max rg 

{QJ.vh 



s.t. (6), 

(12a) 


tr(Qfc) < yPs \\hkf , Vfc, 

(12b) 


Qk p Vfc. 

(12c) 

Next, we 

define a new function F{{Qj.},W) as 



F{{Qk},W)^\^. 

1 +7e 

(13) 

It can be easily shown that the optimal solution {Qk}, W* 
to (PI), is also optimal for (PI') given by 

(PI') 

; max PF) s.t. (12a) — (12c). 

{Q,},W 

(14) 


Hence, we focus on solving problem (Pl^) in the rest of 
the paper. However, since (Pl^ is in general a non-convex 
problem that is hard to solve, we will reformulate it into a 
two-stage optimization problem. Eirst, we constrain the SINR 
at Eve to be %, it thus follows from (13) that F{{Qk},W) 
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is maximized when 7 t, is maximized, which can be obtained 
by solving the following problem: 


(PlM) 

S.t. 


max 

{Qk,m 


PAhoWhol^ 


Ps\gIWho\^ 


cT^glWW^gl + Y.t=i alQkdl 


K 


■=7e, 


(12a) - (12c). (15) 


Let H{je) denote the optimal value of (P1^1) given je. Then 
(Pi') can be equivalently solved by 


(Pl'.2) 


max 

7e>0 


l + g(>) 

1 + 7e 


(16) 


Lemma 3.1: Problem (Pi') has the same optimal value as 
(Pi'.2), and the same optimal solution as (Pl'.l) when % 
takes the optimal solution for (Pi'.2). 

Proof: The proof follows from [7, Lemmas 4.1-4.2}. ■ 

Therefore, (Pi') can be solved in the following two steps. 
First, given any fe > 0, we solve (Pl'.l) to attain 
then we solve (Pi'.2) to obtain the optimal 7*. 


B. Optimal Solution to (Pl'.l) 

Here, we consider solving problem (Pl'.l) by jointly op¬ 
timizing the covariance matrix for the AN at each of the HJ 
helper, Qfs, and the beamforming matrix, W. To facilitate 
the analysis in the sequel, we rewrite the following equations 
in line with our definition of vec(-) [32, Chapter 13]: 


\hlwh:,\^ = 

\vec^ {hohQ)vec{W)\‘^, 

(17) 

h^WW^hl = 

\\hl (g)/vec(W)|p, 

(18) 

\g^Who\^ = 

|vec'^(flro/t-o)vec(W)p, 

(19) 

glWW^gl = 

||g^®Jvec(TY)f. 

(20) 


In addition, tT(W{PshQhQ + ct^J)W^) = ||$u?|p, where 
$ = (/(g)©^)i/2 with 0 = Pshoh^ Hence, (Pl'.l) 

can be rewritten as 


(Pl'.l-RW) : 

max 

W.{QO 

s.t. 


_ Ps\flw\^ _ 

cr2 IIFiruf -f J2k=i hlQkhl + cr^ 

_ Ps\f^w\^ _ 

(T^\\Y2wf+ J2f^^glQ^gl + a‘^ 

(21a) 

(21b) 

(12b), (12c), 


in which w = vec(Ty), fi = vec(h.o^o )’ /2 = 

— T 

YI — Hq (g) J and Y 2 = g'^ ® I. 

As problem (Pl'.l-RW) is non-convex, we define X = 
ww^, Fi F 2 ^ flfl, Fi 4 yf Y 2 ^ 

Y 2 Y 2 and $ = Then by ignoring the rank-one 


constraint on X, (Pl'.l-RW) is modified as 
(Pl'.l-RW-SDR-Eqv) : 


max 

x.{QO 

s.t. 


Pstr{F,X) 


a^Y.X) + 
Pstv{F2X) 


K 

= ie ( Cr^r(y2X) P'^glQkOi 

\ k=l 

tr(¥X) < Pr, 

tr(Qfc) <riPs\\hkf , Vfc, 

^ ^ 0, Qfc A 0, Mk. 


(22a) 

(22b) 

(22c) 

(22d) 


Problem (Pl'.l-RW-SDR-Eqv), via Charnes-Cooper trans¬ 
formation [33], can be equivalently recast as 


(Pl'.l-RW-SDR) : max 

X,{Qk},r 


S.t. cr^ 


Psir{F^X) 

ir{Y,X) + Y. Qk^k + — 1> 

fc=l 

Pst^{F2X) 

= 7e (crltT{Y2X) + J29lQk9l + raA , 


K 




tr(^X) < rPr, 

tr {Qk) < TgPs \\hk\\^ , V/c, 

XPO.Q^PO, V/c, T > 0. 


(23a) 


(23b) 

(23c) 

(23d) 

(23e) 


Lemma 3.2: The constraints in (23a) and (23b) can be 


replaced by (T^tr(yiX) + J2k=i ^kQk^ 

PsXv{F 2X) < fe{(r'^Xv(Y 2 X) + Y.k=i9lQk9l + raf), 
respectively, where both inequalities will be activated when 
problem (PI') obtains its optimum value. 

Proof: See [34, Appendix A]. ■ 

Since problem (Pl'.l-RW-SDR) is a standard convex op¬ 
timization problem and satisfies the Slater’s condition, its gap 
with its dual problem is zero [35]. Now, let A denote the 
dual variable associated with the equality constraint in (23a), 
a associated with the other equality constraint in (23b), /3o 
associated with the transmit power constraint for the AF relay 
in (23c), {/3fc} associated with the transmit power constraints 
for each in (23d), and C, associated with t. Then the 
Lagrangian of problem (Pl'.l-RW-SDR) is given by 

K 

L{Ct) = tr(AJA) + tr(BfcQj,)-f Cr-f A, (24) 

fc=i 


T(T^ < 1 and 


where O denotes the set of all primal and dual variables. 


A = RFi - Aaf Yi - aP.,F 


Bf, = —Xhuh 


k"'k 


C = -Aag 




_ _ Q;7eCr^Y2-, 

afeOlgl - hi, v/c, 

K 

hPr P^'qPshWhkW' 

k=l 




(25) 

(26) 

(27) 


Proposition 3.1: The optimal solution, (X*, {Q^}, t*), to 
(Pl'.l-RW-SDR) satisfies the following conditions: 
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1) rank(Qfc) | ^ 

2) X* can be expressed as 


if/3fc =0, 

if/3^ >0, 


Vfc; 




N^-r, 

n=l 


(28) 


where a„ > 0 Vn, b > 0, rc = rank(C'*) (c.f. (88)) and 

^ G C^* is a vector orthogonal to H = {'nn}nLi 
which consists of orthonormal basis for null(C'*); 

3) According to (28), if rank(X*) > 1, then we have the 
following sufficient condition to yield an optimal solution 
of X with rank-one: 


i:* = (29) 

Qk = Ql Vfc, (30) 

t*=t*+At, (31) 


is also optimal to problem (Pl'.l-RW-SDR), if there 
exists At > 0 such that 


'N?-ra 

X “"t’' 

n—1 





< At < 





n—1 


a„tr (ry^yiJ7„) . 


(32) 


Proof: See Appendix A. ■ 

Note from Proposition 3.1 that if rank(X*) = 1, then 
the optimal w* to (Pl'.l-RW) can be found directly from 
the eigenvalue decomposition (EVD) of X , where X = 
X*jT*. Namely, the upper-bound optimum value obtained 
by solving (Pl'.l-RW-SDR) is tight in this case; otherwise, 
(X*, {Q^},t*) only serves as an upper-bound solution. 

Now, we show that this upper-bound is always achievable 
by a rank-one X*. When rank(X*) > 1, firstly, we check 
whether the sufficient condition proposed in (32) is satisfied. If 
it is met, then a direct reconstruction of (X , {Qf^}, f *) with 

^ 5(t 

rank(X ) = 1 follows according to (29)-(31); otherwise, as¬ 
sume that any optimal solution to problem (Pl'.l-RW-SDR) 
has no zero component, i.e., (X* f 0, {Q^ f 0}, t* ^ 0). In 
addition, the number of optimization variables and the number 
of shaping constraints are denoted by L and M, respectively. 
Since L = X -f 2 and M = X -f 3 for (Pl'.l-RW-SDR), we 
have M < L + 2 satisfied. Thus, according to [36, Proposition 
5.5], (Pl'.l-RW-SDR) has an optimal solution of X that is 
rank-one. Also, the detailed rank reduction procedure based on 
an arbitrary-rank solution has been given in [36, Algorithm 1]. 
Algorithm I for solving (PI') is shown in Table I. 


C. Suboptimal Solutions to (Pl'.l) 

1) Optimal Solution Structure based Scheme: We propose 
a relay beamforming design for (Pl'.l) based on the optimal 
structure of W [37, Theorem 3.1]. First, define iTi = [/tq 9o] 
and H 2 = [ho gg]- Then express the truncated singular-value 
decomposition (SVD) of iTi and H 2 , respectively, as 

H, = XiSiVf, (33) 

H2 = U2'^2V!^. 


TABLE I 

Algorithm I for (PI') 


• Initialize 7e_search = 0 : o : 7 e max and i = 0 

• Repeat 

1) Set i = i -I- 1; 

2) Given 7e = 7e_search(*). 

solve (Pl'.l-RW-SDR) and obtain JT( 7 ^''). 

• Until i = L, where L = [ J + 1 is the length of %_search 

• Set 7 j = 7e_search (argmax| ^ for (PI'.2) 

. Given 7 ^, solve (Pl'.l-RW-SDR) to “obtain (X*, {QI},t*) 
if rank(X*) = 1, apply EVD on X* such that X* = w*w*^\ 
else if the sufficient condition in (32) is satisfied, 
construct (X*, {qJJ, f*) following (29)-(31) and set w* = 
else construct X using the procedure in [36, Algorithm 1]. 
end 
end 

• Recover W* =vec~^{w*) 


Lemma 3.3: The optimal relay beamforming matrix W for 
problem (Pl'.l) is of the form: 

W = uIbU^+ u\c{ui)^, (35) 

where B G and C G unknown 

matrices, and [/f, C/f G satisfy Xf(C/f)-^ = 

I - D'lD'f, C/f (C/f )^ = I - U 2 U 2 , respectively. 

Proof: See Appendix B. _ ■ 

Denote U^ho, U^Hq, U^Qq by ho, ho, Qo, respectively. 
We thus simplify j/igW/rop and \hoBho\'^ 

and [qq Bho\‘^, respectively. Since C has 2{Nt — 2) complex 
variables, we devise a suboptimal design for C to reduce the 
size of variables by_(Xt — 2). Specifically, let C = u'^v^, 
where u' = /rJ/||/ro|| such that = I — u'u'^. 

Hence, ho goWW^gl and (5) can be reduced 

to WB^hlV, WB^SlW^ + [v^u'^^gll^ and Ps\\BhoV + 
a^ti{B^Bf + (^[[vW"^, respectively. Then define b = vec{B), 
= vec{hoho), S 2 = vec(go/io), Y'^ = /ij ® J, = 
g^ ® I, and $' = (/ 0 0'^)^/^ with 0' = Pghoh^ + cr^I', 
Z = bb^, V = vv^, Fi = flf^, F 2 = fill, 
y\ = yI = YfY 2 , and The 

suboptimal design for problem (Pl'.l) by ignoring the rank 
constraints on Z and V is thus given by 

(Pl'.l-subl-SDR) : max Pstr(Fi.Z) 

Z,V,{Q,}.r 

K 

s.t. a-ltT:{YlZ) + ^hlQt.h\FTal = 1, (36a) 

k^l 

Psir{F2Z) < % 

(tr(y2Z) -f \g'^u'^\hY{V)^ + Xfl'iQfefl'I + PO-e j . 

(36b) 

Z) + i7^tr(Z) + a'^tr{V) < rPr, (36c) 

tr (Qk) < TijPs \\hkf , Vfc, (36d) 

T>0,QkhO,yk, ZhO,Vh 0. (36e) 


(33) 

(34) 


Remark 3.1: The variables in (Pl'.l-subl-SDR), i.e., Z G 
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(j^4x4^ V G (j^(JVt- 2 )x(Aft- 2 )^ Qf much reduced size. Fur¬ 
ther, the reconstruction of rt* from V can be briefly explained 
as follows. Given the Lagrangian of (Pl'.l-subl-SDR), the 
KKT conditions with respect to (w.r.t.) V* are given by 

- P*al)I + U*=Q, (37) 

U*V* = 0. (38) 


Post-multiplying (37) with V*, we have (a* 7 e|g'Qtt 




PqCj'^)V* = 0. As a result, if ^ ^ - 




V* = 0; 


otherwise V* = v*v*^, with v* = ^/tT{V*)vo, where 
Vq G c(^t- 2 )xi jg arbitrary vector with unit norm. With V 
solved, (Pl'.l-subl-SDR) reduces to a problem with similar 
structure as (Pl'.l-RW-SDR), and the proof for existence of 
a rank-one Z can be referred to Proposition 3.1. 


2) Zero-forcing: We propose a low-complexity ZF scheme 
for (Pl'.l), in which the jamming signal places a null at 
Bob, and then a semi-closed form solution for W is derived. 
In line with the principle of ZF jamming [17], the jamming 

signal rik is designed as = V k^k such that I — 


H 


hlK 

\\hk\\- 


VkVk , and fik G is an arbitrary random vector, 

fik ~ CA/’(0, Qfe), k = 1,..., AT. Thus, given any W, n^’s 


can be optimized to maximize the effect of jamming at Eve 

by TaaxJ2k=i9lVkQkVk ol, which gives qI = Ckdlol, 
Qk 

~ T 

where and ^ = ^/riP^\\hk\\/J^gk\\ is determined 

by (23d), Vfc. As such, ^^=1 dk^kQk^k 9k turns out to be 
J2k=i 7^!s||/rfc|P||Sfc|P, which is denoted by q. 

With fixed q, (Pl'.l-RW-SDR) can be recast as 


(Pl'.l-sub2-SDR) : max P,tr(FiX) 

s.t. CT^tr(Y iX) -G Ta1 = 1, (39a) 

P,tr(F2X) < % (a^r(F2X) + rg + ral) , (39b) 

tr(¥X) < TPr, (39c) 

X A 0, T > 0. (39d) 


Proposition 3.2: (Pl'.l-sub2-SDR) must yield a rank-one 
solution, i.e., X* = ww*, such that w* = gL>ma,x{Z*), and 

Z* =PsFi - X*a^Yi - a*PsF2 

+ a*fealY2 - /3o*$, (40) 

where i>'max(-) represents the eigenvector corresponding to 
th e largest eigenvalue of the associated matrix, and p = 

, / — (i, . Also, A*, a* and are the opti- 

Y tr($)IZ„,,,(Z*)«zH„(Z*) ^0 

mal dual variables associated with (39a)-(39c), respectively. 

Proof: See Appendix C. ■ 

The only problem in Proposition 3.2 is the dual problem of 

(Pl'.l-sub2-SDR), which admits a much simpler structure to 

solve than the primal one. 


IV. Joint AN-AF Beamforming with Imperfect CSI 

A. Problem Formulation for Imperfect CSI 

We use a deterministic spherical model [25,26] to charac¬ 
terize the resulting CSIs’ uncertainties such that 

Go ={9o\9o =90 + ^9oi^9o^o^9o < !}> (41a) 

Gk ={9k\9k =9k + ^9k, Agf WfcAc/fc < 1}, Vfc, (41b) 
Ho ={ho\ho = ho-\- Aho, Ah^W'^Aho < 1}, (41c) 

Hk ={hk\hk = hk + Ahk, Ah^W'lAhk < l},Vfc, (41d) 
Hk ={hk\hk = hk + Ahk,Ah^W'kAhk < l},Vfc, (41e) 


where gg, gfs, ho, hfs and hfs are the estimates of 
the corresponding channels; Agrg, Agfs, Aho, Ahfs and 
Ahfs are their respective channel errors; the matrices Wo, 
Wk’&, W'q, Wk’s and W'fs. determine the shape of each 
error region. Without loss of generality (w.l.o.g.), we set 
TFo = I/eo, TV' = I/e'^, Wk = I/ek, W'k = 7/4 
and W'l = I/Ck for simplicity, where eg. eg’ ^k, e/, and 
t'l represent the respective size of the bounded error regions, 
k = l,...,K. 

Accordingly, we denote the robust counterpart for (PI') as 


(P2') : 

max min F({Q^},W") 



BoeGo,Bfc65fe,Vfc 


S.t. 

tr (tv (^Pshoh^ + a"/) TV^) < P„ 

(42a) 


tr(Qfc)<i77]s min ||h.fc||^,Vfc, 

(42b) 


hkGHkyk" 



Qk + 0 , Vft. 

(42c) 


An equivalent robust reformulation of (PI'.2) is given by 


(P2'.2) : max 

7e>0 


1 + g(%) 

1 + Fife) ’ 


(43) 


where F{'ye) = 7e and H{je) denotes the optimal value of 
problem (P2'.l) that is given by 


(P2'.l) : 


max _ min 

hQG'H.Q 


s.t. max — 

BfceCfc.vfc 2 

boECo r 


_ P.tr(FiX) _ 

cr/tr(ViX) -f fnQkhk + c'b 

PstT{F2X) _ ^ _ 

tr(Y2^) + Ef=i9fcQfc9l + '^e 


ei 



(44a) 

tr($X) < Pr, 

(44b) 

tr {Q^)<riPs nun ||Afef,Vfc, 

hkGKk.yk 

(44c) 

rank(X) = 1, 

(44d) 

X PO, Qi^PO, Vfc. 

(44e) 


As stated in Lemma 3.1, similarly, (P2') can be proved to 
have the same optimal value as (P2'.2) and the same optimal 
solution as (P2'.l) when 7 e takes its optimal value. As a 
result, (P2') can be solved in a two-stage fashion as well. 
Specifically, given any %, we first solve (P2'.l) to obtain 
H{fe) and then search for the optimal fe to (P2'.2). 
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B. Solutions to (P2M) 

By ignoring (44d), (P2'.l) is recast as 


First, by rearranging terms, (46a) can be equivalently trans¬ 
formed into the following linear form: 


(P2M-RW-SDR-Eqv) : 


max _ min 


Pstr{FiX) 


a^tT{YiX) 
s.t.(44a) - (44c), (44e). 


Ef=i KQkK 




h 

(45a) 


_ mjn PstiiFiX) - 6a^tr{YiX) 

,Vfe 

K 

- 5(jI>Q. (51) 

fc=i 


It is worth noting that due to the rank-one re¬ 
laxation of (P2'.l-RW-SDR-Eqv), solution provided by 
(P2'.l-RW-SDR-Eqv) in general yields an upper-bound for 
which may not be achievable. However, in the sequel 
we insist on solving (P2M-RW-SDR-Eqv) that is regarded as 
an upper-bound benchmark for our proposed problem detailed 
later in this subsection. 

1) 5olMh'ons to (P2'.l-RW-SDR-Eqv).- To make the “max- 
min” objective function of (45) tractable, we first rewrite (45) 
by the equivalent epigraph formulation as 


(P2'.l-RW-SDR-Eqv) : 

max (5 

x,{Qk},s 

s.t. min -—- 

'“So" + Eti K QuK 

(44a) - (44c), (44e). 


P,tr(FiX) 


><5, 

(46a) 

(46b) 


As there are potentially infinite number of constraints in (46a), 
(44a), and (44c), they are semi-indefinite and thus intractable. 
In the following, we equivalently transform these constraints 
to tractable ones using S-Procedure and a generalized S- 
Procedure given in Lemmas 4.1 and 4.2, respectively. 

Lemma 4.1: (S-Procedure [35]) Let oi = 1,2 be 

defined as 


Recalling the following matrix equalities in line with our 


definition of vec(-) operation: 

tr(AB^) = vec^(A)vec(J5), (52) 

vec{AXB) = (A (g) B^)vec(X), (53) 

{A(^B f = (^B^, (54) 

it follows that 

tr{FiX) = h^{ho(^I)X{h^ (^r)h\ (55) 

tr(FiX) = h^{I (g) X)h\ (56) 


3 “ T 

where h G C^‘ = vec(/iQ (g) J). The equivalent channel 

model for h is given hy h = h + t\h, where || Afi f 
(c.f. (41)). By introducing X" = (/iq (g) I)X{hQ ® I) and 
X' = J (g) X, (51) can thus be recast as 

min AfJ'{PsX"-6a‘^X')Ah‘+2U{AfJ'{PsX"-Sa‘^X')h'^} 

K 

- 5(jI>Q. (57) 


fm{x)=X^AmX + 2U{b^x} + Cm, (47) 

where Am = A^ G bm G and Cm G ®., and 

5ft gives the real part of the input entity. Then the implication 
fi{x) > 0 ^ f 2 (x) > 0 holds if and only if there exists 
5 > 0 such that 



(48) 


provided there exists a point x such that fm{x) > 0, m = 1, 2. 

Lemma 4.2: ([38, Theorem 5.5]) The robust block quadratic 
matrix inequality (QMI), 


H 

{F + GX)^ 


C -t- X^B -f B^X + X^AX\ - ’ 

forallJ-X^DX A 0, 
(49) 


Hence, according to Lemma 4.1, the implication ||Ah,|| " < 
AjCg (57) holds if and only if there exists > 0 such 
that the following LMI holds: 


Hi 

Ff 



(58) 


where Hi = PsX" - 6a^X' + w^°'>I, Fi^ = 
{PsX" - 5alX')h‘' and a = h^(PsX" - Sa^X')h^ - 

5Y^^=i^Qk^k ~ ~ Now, (46a) has been 

equivalently reformulated as (58). To further cope with channel 
uncertainties with regards to /r^’s such that (58) holds for 
hk G Hk, k = 1,..., K, we need the following proposition. 

Proposition 4.1: The semi-indefinite constraint of (57) can 
be equivalently recast as the following block matrix inequality: 


is equivalent to 

' H 

F 

G ■ 


0 

0 

0 ■ 


'Hi^^ 

jr(x)H 

f(x) 

ryWH 


'o 

0 

0 

1 

0 ' 
0 

> 0, such that 

F^ 

G 

B^ 

- t 

0 

I 

0 

t 0. 

r^(K)H 

b(^) 

a(k) 

0 

0 

-I 



B 

A 


0 

0 

-D 


.'^1 


^1 J 


L 




(50) 


(59) 



(K) j.(K) 


and are recursively given by 




(fc-i) 


L»(^ 1) |r l)i^ 


Gi' 


(fc-i) 


PsX" -6 (jIX' + I, 


(fc-i)' 


jjik—l) 


k > 1; 
k = 1, 


B 


fc > 1; 

{PsX" -Sa^X')h\ k = l, 


(60) 

( 61 ) 


where H 


ffW = 


_p(fc) ^ 


cf) = h^{P,X" - 6alX')h^ - - 

^llf=k+ifh Q^h\ - - Ef=i^ k = 

1,...,K. In addition, G C(w3+(fc-i)iVt)xJVt ^ 

^ -6Qkhl, Af'> = -6Qk, k = l,...,iG, and 
{^(fc) > 0} denote pertinent auxiliary variables. 

Proof: See Appendix D. ■ 

Next, (44a) is rewritten as 


K 


max 

SkCSk.Vk- 

90 ^Qo 


{PsX" - -ieCrlX') - 7e ^ alQkOl 


k=l 


- leoi < 0, (62) 


where g G ^ ^ = vec(gQ 07) and the equivalent imperfect 
channel model is given by g = g + Ag such that ||Ag|| 2 < 
Nteo- 

Proposition 4.2: The semi-indefinite constraint of (62) is 
satisfied if and only if there exists > 0, fc = 
such that the following block matrix inequality holds: 



1_ 

p(K) 

2 

1- 


piK)H 


p(K)H 



bY 

l(K) 

^2 J 

where 77 2 ^ ( 

fY 

and cY 



0 

0 

0 ■ 


0 

1 

0 

0 

0 

7 

EK -1 


0 , 


(63) 




Fi^'> = 


. (fc —1) .^(fe-1) j 
2 €k-i 2 

^(fc—1) rr(fc —1) 

-PsX" + %a^X' + v(°h, 


fc > 1; 
fc = 1 , 


rj(fe—1) 
jp{k—l) 


fc > 1; 


(64) 

(65) 


i-PsX" + %<J^X')g\ fc = l, 


= g^i-PsX" + %alX')g^ Qj9]+ 

i=i 

7e alQigl + iecrl - v^°'>Nteo - (66) 

z—fc+1 1 — 1 

fc = Also, = feQkOh and 

=%Q^,k = l,...,K. 


Proof: See Appendix E. ■ 

Last, we rewrite (44c) to facilitate the robust optimization 
against the errors introduced by Ahfs. By applying Lemma 
4.1, (44c) holds if and only if there exists Pk ^ 0, k = 
1,... ,K, such that the following LMI constraint is met: 

gPsI + Pkl gPshk ^ ^ 

. H - AO, Vfc. 

tiPsK gPs\\hk\\i-tT{Qk)-Pk^'k 

(67) 

As such, (P2'.l-RW-SDR-Eqv) is now simplified as 
(P2'.l-RW-SDR-Eqv) : max 5 

X,{Q,}.5 

s.t. (59), (63), (67), (44b), (44e). 


Because of the non-convex term such as 5X' in (59), 
problem (P2'.l-RW-SDR-Eqv) remains very hard to solve. 
We thus use the bisection method [35] w.r.t. 6 to solve 
it. However, using bisection in addition to solving (P2'.2) 
by one-dimension search over fe may lead to very high 
complexity. As a result, we propose an alternative problem 
to approximate H{fe)- 

2) Solutions to (P2M-RW-SDR).- We propose to approx¬ 
imate H{fe) by the optimum value of the following problem. 


(P2M-RW-SDR) : max min P^tr(FiX) (68a) 

ho gPo 


K 


'tTiYiX) + YhkQkK+Tal<l, (68b) 

k=l 

Pstr{F2X) 


s.t. _ max 

h^<=Hi^,Vk 
h-Q&'UQ 

max - 

'^MY2X) + ELi alQkol + ■ 

tr(¥x) < TPr, 

ir{Q,,)<TgPs min \\hk\f, ^k, 

hk^HkS/k 

^ A 0, Qfc A 0, Vfc, T > 0. 


< 7e. (68c) 

(68d) 

(68e) 

(68f) 


Remark 4.1: It is worth noting that as the numerator and the 
denominator of the objective function in (P2'.l) are coupled 
by common uncertainty Hq, Charnes-Cooper transformation, 
in general, cannot be applied to realize equivalent decoupling. 
As a result, (P2'.1-RW-SDR) yields a more conservative ap¬ 
proximation for H{fe) than (P2'.l-RW-SDR-Eqv). However, 
considering that (P2'.1-RW-SDR) needs to be solved only 
once for given % in contrast with (P2'.l-RW-SDR-Eqv) 
requring isection over 5, we exploit it in the sequel. The 
effectiveness of this approximation will be evaluated in Sec¬ 
tion V-B. 

To proceed, we rewrite (P2'.1-RW-SDR) as 

(P2'.1-RW-SDR) : max 5 

X,{Qj,}.5.r 

s.t. min Pstr(ifiX) > <5, (69a) 

ho^i-lQ 

(68b)-(68f). (69b) 

Lirst, by rewriting F = f [fi , where fi=fi + A/^, in 
line with Lemma 4.1, the implication jjA/^lP < ||h.o|peo ^ 
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(69a) holds if and only if there exists s'(0) > 

0 such that the 

following LMI constraint is satisfied; 


PsX + s'(o) J PsXfl 


0 . 


_ (70) 

Next, as tr(l^iX) = y'{X'y\ (c.f. (56)), where y^ = 
vec{hQ (SiI), after some manipulation, (68b) holds if and only 
if there exists > 0 such that 


s"(0)j_a2X' -a^X‘ 
[ -c^ly^X' c 


Vi 


^0, 


(71) 


where c = -aly'{X'y\ - Y.k=i - t(jI + I - 

s"(^)Nte'Q. Then (68b) can be rewritten as 

(71) for hk G Tffc, 'ik, (72) 

which is handled by the following proposition. 

Proposition 4.3: The semi-indefinite constraints in (72) can 
be replaced by the following LMI constraint; 


0 , 


rfjiA 

p(K) 

^iK)( 


0 

0 

0 ' 

p{K)H 

cW 

p(K)H 

- .f'A) 

0 

1 

0 

q(K)H 

piK) 



0 

0 

-I 


where and are recursively given by 


(73) 






^(fc-1) jjy 

k = 2,...,K; 

[fjW = = -a^,X'y\, 


r(fc-l) 




■^ik-iy 

pik-l) 


where = G‘i^\ = -Q^hl, A''^’ = -Q^, = 


,(fe) 


t 


(74) 

A) = 


■t 


-atVi X'yl - Q^h< - J^tk+ihl QM - raf + 1 - 

fc = 1,..., iT, and > 0} denote 

the auxiliary variables. 

Proof: See Appendix F. ■ 

Proposition 4.4: The constraint in (68c) is guaranteed if 
and only if there exists > 0, k = 1,... ,76, such that 
the following LMI holds; 


^0, 

(75) 


vH^k) 

fA) 

G(k) - 


0 

0 

0 ■ 

p(K)H 

c'iK) 

pdK)H 


0 

1 

0 

q(K)H 

B'A) 

A'0<) _ 


0 

0 

-I 

€.K -1 


where and are recursively given by 

r A'(k-i ), 

H^k)^ 


r y(fc-i) 1 s'"-!’/ 

~ efc_i 

Q(k-1)H' 


pf{k-l) 

Q{k-1) 

Jjik-1) 

1 -*7 — 

p{k-l) 


k = 2,...,K; 

= -PsX" 
= {-PsX'‘ 


- feolx' 


s(0)l, 


■leO, 


■X') g\ 


(76) 


in which = %Qk9\i = leQk^ 


„i(k) — 


= g^{-PsX" -f + 7 e Ei=l sJQjSj 


le'}lf=k+l 9 jQ^ 9 \ + leTdl - S^^Ateo - ^ = 

1,..., 76, and > 0} denote the auxiliary variables. 

Proof: It is observed that (68c) differs from (44a) in the 
only respect that is replaced by rCg. Hence the proof for 
Proposition 4.2 can be directly applied herein by substituting 
rcTg for CTg. ■ 

Last, by replacing “yPs” in (44c) with “ryPs" in (68e), 
(68e) can be replaced by a similar LMI as (67), denoted by 
(68e'), in which the pertinent auxiliary variables are denoted 
by {ptk > 0}. 

Consequently, the equivalent reformulation for problem 
(P2'.1-RW-SDR) can be summarized as 


(P2'.1-RW-SDR) : max 5 

70 ),„'( 0 ),„"( 0 )_ 

{^fc} 

s.t. (70), (73), (75), (68e'), (68d), (68f), 

> 0, s'(°) > 0, s"(°) > 0, (77a) 

g(k) > 0, > 0, /Xfc > 0, Vfc. (77b) 


C. Proposed Rank-One Solutions to (P2'.l) 

(P2'.1-RW-SDR) is convex and can be solved efficiently by 
convex optimization tools such as CVX. Next, we derive the 
Lagrangian of (P2'.1-RW-SDR). Note that in the following 
expression, we only consider the uncertainties regarding Hq, 
hk’s, hfs, Qq and gfs when 76 = 1 for the purpose of 
simplicity and the results can be easily extended to the case 
of 76 > 1. Denote the dual variables associated with (68d), 
(70), (73) and (75) by (3q, W, V and Y, respectively. Then 
the partial Lagrangian of (P2'.1-RW-SDR) w.r.t. X is 

L{n) = tr(AX), (78) 


where FJ is the set of all primal and dual variables, and 

A = PsWi,i + 2P,JR{/VfE + PsW2,2Fi 


Nt 

X(v<:;>+2»{v;“>}+F<« 


2=1 


- PsK ® I)yiAho <S) I) + 2aAe ^ } (79) 

2=1 

-2Ps^{{h^ (Sl)Y2,i{ho (Si)} 


Nt _ 

- Psy 2 , 2 iho ® 1 ) 9 ^ 9 ^(ho ® /) + J 2 ^S 2 - 


In (79), Fi = f\fF, Wt,j, i,j = 1,2, Vt,j, z,j = 1,... ,3 
and Yij, i,j = 1,...,3 are the block submatrices of 
g £;(JV 2 -l-l)x(Arf-|-l)^ y g £{N^+Nt-\-l)x{N^+Nt-\-l) 

Y G (C,AA‘^Nt+i)x{Nf-i- 2 Nt-i-i) jjjg same size as block 
submatrices in (70), (73) and (75), respectively. Moreover, in 
(79), V'i,2 = y\Vi2t y2,2 = V2,2y\yi, Y2,1 = AyT,2 
and Y 2,2 = 2 / 2 , 29 ^ 9 ^- Furthermore, , y^i '2 ^ 2 * 2 ^ 

are the zth block diagonal submatrices of G 
1^1,2 G cA^n} y 2_2 G respectively; 
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and Y 2 2 are the zth block diagonal submatrices of 1 ^ 2,1 G 
and y^ 2,2 G respectively. 

Proposition 4.5: 1) The optimal X* to 

(P2'.1-RW-SDR) is expressed as 

X*= Y. anfinfin +blt, (80) 

n—1 

where d„ > 0, Vn, 5 > 0, and | G ^ ^ is a unit-norm 
vector orthogonal to H (c.f. (28)). 

2) According to (80), if rank(X*) > 1, i.e., there exists at 
least one d„ > 0, we reconstruct a solution to problem 
(P2'.1-RW-SDR) using 

X* = (81) 

5* = (5*, (82) 

T*=T*, (83) 

^ 5(5 

while {Qf^} are obtained by solving the following feasi¬ 
bility problem provided that X , 6*, and r* are given by 
(81), (82) and (83), respectively: 

(P2'.l-RW-SDR-sub) : max 0 

{Qfc},s"(0), 

{<!"('“)},{Mfc} 

s.t. (73) given X ,f*, (68e') given f*, 

Qk ^ 0, /Xfc > 0, Vfc, 

> 0, > 0, Vfc. 

Proof: See Appendix G. ■ 

The scheme that solves (P2') is summarized in Table II. 

TABLE II 

Algorithm II for (P2') 

. Initialize 7e_search = 0 : q' : and i = 0 

• Repeat 

1) Set i = i + l\ 

2 ) Given = 

solve (P2'.1-RW-SDR) and obtain 

• Until i = L', where L' = + 1 is the length of 7 ' 

. Set 7 ; = 7'_,earch ^argmax^ }) (P2'-2) 

. Given 7 J, solve (P2L1-RW-SDR) to'obtain (X*, {QJ}, <5*, t*) 
if rank(X*) = 1, apply EVD on X* such that X* = w*w*^\ 
else 

- construct {X*,&*,t*), according to (81)-(83) and set w* = 
s/bl', 

— given X , S* and r*, 

obtain {Q*f\ by solving (P2'.l-RW-SDR-sub). 

end 

• Recover W* =wec~^(w*) 


V. Numerical Results 

Here we provide numerical examples to validate our results. 
We assume a typical scenario where the K helpers are evenly 
distributed around Alice with a radius of pk = 2m and 9k = 
(radian by default), where 9k is the angle of direction 
(w.r.t. the Alice-relay link by default) of the fcth helper, k = 


1,... ,X. Alice, Bob and Eve are, w.l.o.g., assumed to have 
the same distance away from the AF relay with their angle 
of direction tt, tt/O and IItt/O, respectively. We also assume 
channel models with both large-scale fading, i.e., path loss and 
shadowing, and small-scale fading, i.e., multi-path fading. The 
simplified large-scale fading model is given by [39] 

( 7 \ —a 

— j , for d > do, (84) 

where z is a log-normal random variable capturing the effect 
of shadowing with the standard derivation cr = 4dB, Aq = 
10“^, d is the distance, do is a reference distance set to be 
Im, and a = 2 is the path loss exponent. Specifically, the 
channels including fife’s, fig, fig and g^, are assumed to suffer 
from Rician fading while the channels from the HJ helpers to 
Bob (fife’s) and Eve (flffe’s) follow Rayleigh distribution due 
to the missing of line-of-sight (LOS) components with their 
respective average gain specified by (84). Take fife, Vfc, as 
an example, fife = -^^^fife + where fife is 

the LOS component with ||fife ||2 = D (c.f. (84)), fife is the 
Rayleigh fading component denoted by fife ~ CAf {0, DI), 
and Kn is the Rician factor set to be 3. Note that for the 
involved LOS component, we use the far-held uniform linear 
antenna array to model the channels [40]. In addition, unless 
otherwise specihed, the number of HJ helpers, K is set to be 
5; the AF relay is assumed to be 5m away from Alice; the 
EH efficiency, ry = 0.5 and = —50dBm. The 

results presented in Section V-A are obtained by averaging 
over 500 times of independent trials. 

A. The Perfect CSI Case 

We compare the proposed optimal solutions with three 
suboptimal schemes in the case of perfect CSI. One suboptimal 
scheme, denoted by “Suboptimal 1”, is introduced in Sec¬ 
tion III-Cl by exploiting the optimal structure of W. The other 
described in Section III-C2 is known as optimal null-space ZF, 
denoted by “Suboptimal 2”. Specihcally, each jamming beam 
life is restricted to lie in the orthogonal space of fi^ such 
that life’s cause no interference to the IR while maximizing 
its effect of jamming at the eavesdropper. As a benchmark, 
we also present the well-known isotropic jamming that is 
particularly useful when there is no Eve’s CSI known at each 
HJ helper, Hfe, Vfc [41], denoted by “Suboptimal 3”. Note that 
the difference between “Suboptimal 2” and “Suboptimal 3” 
only lies in the design of jamming noise, for which the former 
also aligns the jamming noise to an equivalent Eve’s channel to 
confront Eve with most interference, while the latter transmits 
isotropic jamming with fik ~ CA/'(0, r 7 Ps||fife|p/(Xi — 1)), 
k = 1,..., K, in directions orthogonal to fife’s, due to lack of 
knowledge of Eve’s channel and thus is expected to be less 
efficient than “Suboptimal 2” with perfect CSI. 

First, we study the secrecy rate at the receiver versus the 
transmit power of the transmitter, Pg with Pr = lOdBm. Fig. 2 
demonstrates that for both cases of Xj = 3 and Nt = 5, 
the average secrecy rate increases and tends to be saturated 
as Pg goes to 30dBm. It also illustrates that “suboptimal 1” 
and “suboptimal 2’’ closely approach the optimal solutions 
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Source transmit power(Pg) (Watt) 



Fig. 2. Secrecy rate versus Alice’s transmit power with perfect CSI. 



Relay transmit power(P^) (Watt) 

Fig. 3. Secrecy rate versus the relay’s transmit power with perfect CSI. 

while “Suboptimal 3” is outperformed more succinctly with 
larger number of antennas at the AF relay and the HJ helpers. 
Moreover, with Nt increasing, the average secrecy rate gets 
larger as a result of the higher array gain of the AF relay and 
more available power transferred to the HJ helpers. 

In addition, we show in Fig. 3 the secrecy rate achieved by 
different schemes versus the transmit power of the AF relay, 
Pj. with Pg = 30dBm. It is seen that the average secrecy rate 
first grows faster and then slower, since when Pr increases, 
not only the desired signal but also the noise yielded from 
the first transmission phase is amplified to a larger extent. 
In addition, the performance gap between the optimal scheme 
and suboptimal schemes is almost negligible. Similar to Fig. 2, 
“Suboptimal 3” appears to have certain performance loss from 
the optimality but is considered as a promising scheme when 
no Eve’s CSI is available at the HJ helpers. 

B. The Imperfect CSI Case 

Now, we consider the imperfect CSI case and compare the 
proposed scheme Robust SDR with HJ, which is obtained 
by solving (P2'.l-RW-SDR-sub), against some benchmarks. 
Note that there are two upper-bound benchmark schemes, 
namely. Robust SDR with HJ and Robust-eqv with HJ, as well 
as two lower-bound benchmarks, which are Robust w/o HJ 
and Non-robust with HJ. For Robust SDR with HJ {Robust- 
eqv with HJ), given any %, H{fe) is approximated by 


Fig. 4. CDFs of the achievable secrecy rate. 

solving the rank constraint relaxed problem (P2'.1-RW-SDR) 
((P2'.l-RW-SDR-Eqv)). On the other hand, for Robust w/o 
HJ, we solve (P2'.1-RW-SDR) by setting Q^. = 0, Vfc while 
for Non-robust with HJ, (11) is evaluated by applying the 
optimal solutions to (Pl'.l) assuming perfect CSI, to the 
actual channels including errors that are generated from the 
sets dehned in (41). 

To assess the worst-case secrecy performance, we use the 
metric, namely, secrecy outage probability, defined as [42]: 


p = P^{r < r^), 


(85) 


where rj] obtained by solving (P2') is termed as the 100p%- 
secrecy outage rate. 

The parameters are set identical to those in the perfect CSI 
case. Regarding the uncertainty model in (41), we introduce 
the uncertainty ratios associated with eo. Cq’ ^'k 
ao, ttg, ttfe, aj; and respectively. For instance, ao is 


a 


2 _ 
0 ~ 


eo 

nwaoW^r 


( 86 ) 


while ag, a^’s, a^’s and a'l’s are similarly dehned and thus 
omitted here for brevity. Besides, it is reasonable to assume 
that the channel estimates w.r.t Eve suffer from more errors 
than those for Alice and Bob. Hence, we set = aj? = 
= 1% while = al = 10%, k = 1,...,K, unless 
otherwise specihed. 

Eig. 4 demonstrates the cumulative density function (CDE) 
of the achievable secrecy rate from 1000 samples of ran¬ 
dom channel errors uniformly distributed over the sets de¬ 
hned by (41) given hxed actual channel realization. We set 
Pr = 20dBm, Ps = 30dBm, Nt = 3, K = 5 and 
aQ = aj? = = 2%, k = 1,...,K. Despite being 

suboptimal to the upper-bound schemes of “Robust SDR with 
HJ” and “Robust-eqv with HJ”, the proposed “Robust with HJ” 
scheme outperforms its non-robust counterpart “Non-robust 
with HJ” particularly in the low range of probability, and 
overwhelmingly surpasses the “Robust w/o HJ”. Eor example, 
“Robust with HJ” can achieve a secrecy rate of around 
3.5bps/Hz in the 3% worst case versus that of 3.3bps/Hz and 
l.Obps/Hz for the “Non-robust with HJ” and “Robust w/o HJ”, 
respectively. The solutions for “Robust SDR with HJ” is also 
seen to admit very little gap from those for “Robust-eqv with 
HJ”, which suggests that approximating H{fe) by solving the 
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Fig. 5. Secrecy outage probability for K = S and A" = 5 HJ helpers, 
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Fig. 6. Secrecy outage rate versus the normalized channel errors. 

complexity reduced “Robust SDR with HJ ” leads almost no 
performance loss. 

Fig. 5 illustrates the CDF of the achievable secrecy rate from 
1000 samples of random channel errors generated in the same 
way as Fig. 4, with Pj. = 20dBm, Pg = SOdBm and Nt = 3. It 
is observed that proposed solutions to “Robust with HJ” nearly 
achieve their upper-bound rank constraint relaxed solutions, 
i.e., SDR, to “Robust upper SDR with HJ” throughout the 
whole range of outage probability. Moreover, the “Robust 
w/o HJ” yields the worst performance. In particular, when 
the outage probability falls to 3%, the “Robust w/o HJ” 
achieves a worst-case secrecy rate of less than Ibps/Hz while 
the proposed scheme can still guarantee an outage rate of 
rough 1.64bps/Hz and 2.07bps/Hz for K = 5 and K = 6, 
respectively. Also, it is observed that increasing the number 
of HJ helpers will improve the secrecy performance, but we do 
not draw conclusions on the extent to which the secrecy rate 
can increase, since it also depends on the level of channel 
estimation inaccuracy. For example, more HJ helpers may 
also yield larger interference to the legitimate receiver if the 
channels from HJ helpers to Bob are not as well estimated 
as this instance of = 1%, Vfc. Hence we suggest that in 
practice, a mild number of HJ helpers are sufficient in view 
of the trade-off between complexity and performance. 

Fig. 6 shows two different levels (p = 0.20 and p = 0.30) 
of secrecy outage rate versus the channel uncertainty ratios 
(assuming ao = ak, k = 1,..., K), in which P^ = 30dBm, 


Fig. 7. Secrecy outage rate versus the relay’s transmit power. 

Ps = 30dBm, Nt = 3 and AT = 5. It is observed that the 
secrecy outage rate by the proposed schemes decreases slowly 
with the eavesdropper’s CSI error ratios, which validates the 
motivation of the worst-case robust optimization. It is worth 
noting that the advantage of the HJ protocol is more significant 
when the normalized channel uncertainty of Eve’s channels 
surpasses 10%, since the HJ scheme provides more degree 
of freedom for robust design and thus capable of guaranteeing 
larger worst-case secrecy rate against worse channel conditions 
compared to that without HJ. The reasonably suboptimal 
performance of the proposed “Robust with HJ” is also seen as 
from Figs. 4 and 5. 

Fig. 7 studies the 100p%-secrecy outage rate for p = 0.05 
and p = 0.20, respectively, versus the transmit power of the 
AF relay. Specifically, we set Pg = 30dBm, Nt = 3, and K = 
5. As observed similarly from Fig. 6, the robust schemes with 
the assistance of HJ helpers perform considerably better than 
solutions without HJ helpers. Furthermore, when the transmit 
power is set relatively large, i.e., Pg = 30dBm, it is seen that 
continuously increasing P^ does not contribute much to the 
secrecy performance, because in this situation the increased 
amplified noise at the AF relay compromises the performance, 
which provides useful insight for practical setting of Pr. In 
addition, the proposed “Robust with HJ” is observed striking 
a good trade-off between optimality and complexity compared 
with the two upper-bound solutions. 

VI. Conclusion 

This paper considered improving the secret wireless com¬ 
munications in a multi-antenna AF relay wiretap channel via 
a novel harvest-and-jam (HJ) relaying protocol. The AN co- 
variance matrices at HJ helpers and the AF relay beamforming 
matrix have been jointly optimized to maximize the achievable 
secrecy rate and/or worst-case secrecy rate at the legitimate 
receiver subject to the transmit power constraints of the AF 
relay as well as the HJ helpers, on perfect and imperfect CSI 
occasions, respectively, using the technique of semi-definite 
relaxation (SDR). The SDR was shown tight for the perfect 
CSI case while suboptimal rank-one reconstruction algorithms 
for the robust formulation under imperfect CSIs were pre¬ 
sented achieving promising tradeoffs between complexity and 
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performance. The effectiveness of the proposed schemes were 
also verified by numerical results. 

Appendix 

A. Proof of Proposition 3.1 

The KKT conditions of (Pl'.l-RW-SDR) are given by 

A*X* = 0, (87a) 

BlQl = 0yk, (87b) 

P: {tiiQD - T*pPs\\hkf) = 0,Vfc. (87c) 

According to (26), if for certain k, fit = 0, then Bl = 

~ * ~T ^ 

—X*hj.hj. +a*jegk9k rank(B^) < rank(/i^,/i^ ) + 

rank(gr^g^) = 2, which yields rank((5^) > Nt — 2 as 
a result of (87b). Otherwise, when > 0, we will have 

rank(Bfc) > rank(-/3*J - X*hlhl) - rank(a*%g*c/^) = 
A^t — 1 [7, Lemma A.7], which implies rank((5^) < 1. 
However, rank((5^) cannot be 0, since otherwise tr(Q^) — 
PrjPsWhkW^ < 0 and thus /3^ = 0 according to (87c), which 
contradicts to > 0. Hence, when > 0, rank(Q^) = 1. 

Next, define C* = -X*a^Yi - a*PsF 2 + a*%a^Y 2 - 
/3q $ and according to (25), we have 

A* = PsFi + C*. (88) 

Then define r^, H and rj„, n = 1,... (c.f. (28)). 

Similar to the approach used in [7, Appendix B], we discuss 
the structure of the optimal X under two cases. 

(1) Case I: Tc = N^. As C* is full-rank, rank(A*) > rc — 

1 = Nf — 1 and hence — 1 < rank(A.*) < N^. If 
rank(A.*) = Nf — 1, rank(null(A*)) = 1 and it follows 
that X* = by assuming ^ as the only basis of 

null(A*). Otherwise, according to (87a), we obtain X* = 
0, which ceases the secrecy transmission and cannot be 
the optimal solution to (Pl'.l-RW-SDR). 

(2) Case II: rc < Nf. If C* is not full-rank, rank(A*) > 
Tc — 1. Then by pre-multiplying and post-multiplying 

S H with both sides of (88), we have 

Vn ^*rin = PsVn ^iVu + Vn = P-sVn Vn. 

(89) 

According to (24), it is necessary for A* ^ 0 to obtain 
an optimal solution of X* and therefore r 7 ^A*? 7 „ < 0, 
which conforms to PsT]^F iij^ > 0 if and only if 
A* = 0 and = 0. Hence, H C null(A*), i.e., 

— rank(A*) > N"} — Vc ^ rank(A*) < Tc- Next, we 
show rank(A*) Tc by contradiction. If rank(A*) = Tc, 

H = null(A*), and X* = J2n=i^'' ‘^■nVnVn ■ However, 
in this case, since f irj„ = 0, Psti{FiX*) = 0, which 
is apparently not optimal. Hence, we have rank(A.*) = 
Tc — 1 and thus rank(null(A*)) = — rc + 1. This 

indicates that besides the basis in H, null (A.*) spans over 
an extra dimension of basis, which is denoted by and 
hence X* = 

Assume that (X*, {Q^}, t*) is the optimal solution to 
(Pl'.l-RW-SDR) with rank(X*) > 1. Then construct a new 


solution {X , Qfc, t*} according to (29)-(31). Now, we check 
if the reconstructed solution is feasible if (32) holds. First, 

K 


X" ^ „ 

1- ' I ^2 

k+T CF^ 


a^triYiX ) + J 2 KQkh 

k^l 

(_ ( 

< CT^tr I Yi X* - ^ 

\ n^l 

_ 

(^ntr{Yir]„r]^) 


K 


k 


J2~^kQfh 








K 


(a) 




= aMyiX*) + J2^kQlK + T*at < 1 . 
Moreover, 

P,tr(F2X*) = P,tr I F2(X* - ^ 


(90) 


Nf-r 


n—1 


K 


< 7e (altT{Y2X*)+J29kQk9l^r*al 




Nf-r a 


+ Je (TeAT-Cr^tr X 2 


anVnV 


H 


K 


= le\(T, 


'MY2X) + Y.9lQk9l 


I 2 

+ T CT. 


(91) 


fc=i 


In addition, (23c)-(23e) are easily shown to satisfy. In the 
above, (a) and ( 6 ) hold due to the feasibility in (23a) and 
(23b), respectively. Further, Pstr(F’iX ) = Pstr(F’iX*) 
shows that the reconstructed solution achieves the same opti¬ 
mum value as that of (Pl'.l-RW-SDR). Hence, an optimal 
solution to (Pl'.l-RW-SDR) with rank-one X is ensured. 


B. Proof of Lemma 3.3 
First, we construct W as 


W = [U\,{Ui)^ 


B 

D 

t/ 


C 

E 


[U2,U^f 


= U\BUf + U\C{U^f + {UifDU^ 


H 


(92) 


where B G C^x^, C G C2x(^*-2), D G C(^‘-2)x2 and E G 

C(^*- 2 )x(JVt- 2 ) 

are undetermined matrices. Then according to 
(33) and (34), it follows that |/io TF/iqP = JSLTf/ iqP 

and hgWW^hl = \\B^U'^hl\\^ + \\C^U'^hl\\^. Sim¬ 
ilarly, we also have IqqWHoI'^ = \gQU\BUf and 
glWW^gl = \\glU\Br + \\glu\C\\\ Thus, 7 ^ (c.f. (9)) 
and 7 e (c.f. (10)) do not depend on D and E. 

Next, by substituting (92) for W in (5), we have Pr > 
PsiWBU^hoW^ + WDU^hoW^) + a^tT{B^B + C^C + 
D -f E^E). Since (PI') is a secrecy rate maximization 
problem subject to the given P^, it turns out that given the 
optimum secrecy rate, Pr is the minimized required power 
by taking D = 0 and E = 0, while B and C cannot be 
determined directly. Thus, W = U\BUif + U\C{U^)^. 
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C. Proof of Proposition 3.2 

Denoting the dual variable associated with (39a), (39b) 
and (39c) by A, a and /3o, respectively, the Lagrangian of 
(Pl'.l-sub2-SDR) is expressed as 

L{x) = 

tr - aP^F^ + - /3o^)X) 

+ ( —AcTj + afe{<l + CTg) + PoPr) T + A, (93) 

where x = ct)/3o} denotes the set consist¬ 

ing of all the primal and dual variables. Since problem 
(PlM-sub2-SDR) satisfies the Slater condition, its optimum 
value admits zero duality gap with its dual counterpart. Fur¬ 
thermore, according to (93), in order for the dual function to 
be bounded from above, the following constraints must hold; 

Z = PsFi - Acr^Y 1 - aPsF2 + a%(J^Y2 - /3o^ ^ 0, (94) 
- \al + afe{q +crl) + PoPr <^- (95) 


[35, A. 5.5], for (97), we have 

■ffib) Fib) 

Ffm fP-yjW 


Gib) ■ 
[Bib)« 
AD/ 


(^Aib) + t;i^^ ^ [Gib)« Bib)] AO, 


(1) 


which can be reexpressed as 

rAib) + t;i^ Gib)« 
Gib)'" //lb) 

Fib)^^ Fib)^ 


Bib) 

Fib) 


AD 


|,(D 


A 0. 


(98) 


(99) 


Next, assume that the robust design for (58) has been 
considered against the precedent fc — 1 uncertainties, i.e.. 


■//ib-D 

Fib- 1 ) 

Gib-D 1 


0 

0 

0 

Ffk- 1 )H 

Afc-D 

^1 

to 

1 

(fc- 1 ) 

0 

1 

0 

T 


Bib-D 

Aib-D _ 


0 

0 



AO, fc > 2. (100) 


The dual problem is therefore given by 


(D-Pl'.l-sub2-SDR) : min A 

\,a,Po 

s.t. (94),(95), (96a) 

(A,a,/3or>0. (96b) 


It is observed that Z is of the same form as the Hessian 
matrix with respect to X without rank relaxation. Accord¬ 
ing to [43, Theorem 2.1], Z 0 implies that the SDR 
problem (Pl'.l-sub2-SDR) is tight in this case, i.e., 3w* 
such that X* = w*w*^. Moreover, since KKT condition 
necessitates Z*X* = 0, it follows that w* is the eigenvector 
corresponds to the zero-eigenvalue of Z*. Hence, we have 


w* = p.v^^^{Z*), where p = .. ^ „ 

due to the power constraint of (23c), which completes the 
proof. 


D. Proof of Proposition 4.1 

First, given hk, k = 2,... ,K, fixed, only consider the 
uncertainty of hi. Since ||A/ii ||2 < e", we have 1 — 
(Ah-A ^ Ah^ > 0. By applying Lemma 4.2 to (58) with 
= PsX” - SalX' -f rt;(o)/, = {PsX" - 

Sa^X')h\ = 0, c^i^ = h^{PsX" - 5alX')h' - 

ihlQihX - Bi(i) = 

—5Qih\, and Ai^) = —there exists tub) > 0 such 
that the following LMI holds: 


r 

Fib) 

Gib) - 


0 

0 

0 ■ 

Fib)-w 

cp) 

B/Dw 

-lyb) 

0 

1 

0 

T 

_Gi(Dw 

Bib) 

Alb) _ 


0 

0 

—1 


Note that for Qi A 0, there always exists ryb) > Q such 
that + Alb) 0 and we assume that such constraint 

is applied. According to the property of Schur-Complements 


Applying a similar procedure as that for (97), (100) can be 
recast as 


rEE^ + Aib-D 

^fe-1 

Gib-D 


//ib-D 


Bib-D 

Fib-i) 


Afc-D 


„(fc-D 


A 0. 

( 101 ) 


Then given hi, i = k + 1,..., K fixed, accommodate the 
fcth uncertainty, i.e., hk G TLk, for (101). By applying Lemma 
4.2 to the uncertainty of hk, the implication || < ej] 

(101) holds if and only if there exists ryb) > Q such that 


-//lb) 

Fib) 

Gib) 1 


0 

0 

0 - 

Fib)-^ 

cb) 

Bi(Dff 

(fc) 

0 

1 

0 

T 

.Gib)« 

Bib) 

Alb) _ 


0 

0 

-I 

'Ff- 


where 


//lb) = 
Fjb) = 


Aib-D + 


A'^-D 


Gib-D 


Bib-D' 

Fib-i) 


Gib) 


//ib-D 

= 0, 


(103) 


cb) = h^{p,x'' - 6<j^x')h^ - Sj:j=ihjQjh] - 

SElk+ifhQ^hl-6al-w(^^Nie',-EtiwPK Bi^ = 
—5Qkh\ and Aib) = —6Qk, k > 2. Thus, using the 
method of mathematical induction, (58) holds for hk G Hk, 
k = 1,... ,K, if and only if there exists {ty(fc) > 0}, such 
that (59) is satisfied, which completes the proof. 


E. Proof of Proposition 4.2 

Taking the similar procedure as that for dealing with (57), 
the implication ||Ag|p < Ajeo ^ (62) holds if and only if 
there exists yb) > g such that the following LMI holds: 



F2 

C2 


A 0, 


(104) 


where H 2 = -PsX" + %afX' + yb)j, F 2 = 
{-PsX" + %a^X')g^ and C 2 = {-PsX"+ %a^^X')g^ + 
> J 2 k=i 9lQk9l+fe-'^^°^ (44a) has been equivalently 
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reformulated into (104). Then, given fc = 2,..., iT, fixed, 
applying similar procedure to that in Appendix D, it follows 
that there exists > 0 such that the following LMI holds: 


\hP 


G2('> 1 


0 

0 

0 ■ 


cP 

<-2 

_B2(bw 


0 

1 

0 

T 



A2(1) . 


0 

0 

-1 

€1 - 


Since 0 always holds, (105) is equivalent to 


the following LMI: 


A 2 


( 1 ) 




B20)» 




J?2W 

41 ) 


.,( 1 ) 


^ 0 . 


( 106 ) 


Next, devising the method of mathematical induction again as 
that for (99), (104) holds for G Vfc, if and only if there 
exists {v(k) > 0}, such that (63) is satisfied, which completes 
the proof. 


F. Proof of Proposition 4.3 

We only sketch the proof herein since it is quite similar 
to that of Proposition 4.1. First, apply Lemma 4.2 to (71) 
given /Tfc’s, k = 2,... ,K, fixed and obtain an initial LMI. 
Next, manipulate the resulting LMI according to the property 
of Schur-Complements to facilitate using Lemma 4.2. Then, 
repeat this procedure until all the semi-indefinite constraints 
w.r.t. hfs have been incorporated into an equivalent LMI. 


G. Proof of Proposition 4.5 

According to the KKT conditions of (P2'.1-RW-SDR), we 
have A X* — 0 , where A is given by (79). Define C = 
A —W2 2PsFi with rank(C ) denoted by fc. Then take the 
similar procedure as Case I and Case II in Appendix A, it 
can be obtained that X* = ^nVnVn + 

Next, we prove the second half of Proposition 4.5. Accord¬ 
ing to (81), 

PstiiFiX*) = PstiiFiX*) >_min Pstr(FiX) > 6*, 

ho gHo 

(107) 


and thus (69a) holds true, which implies that the same 
optimal value as (P2'.1-RW-SDR), i.e., d*, is achievable. 
However, since the constraint in ( 68 c) is ignored, the global 
optimal 7* for (P2'.2) via solving (P2'.1-RW-SDR) is 
probably violated in (P2'.l-RW-SDR-sub). For example, 
- ^5 -= 70 ^ p(^*'] which results 

in the actual objective value for (P2'.2), smaller than 


l + g(7e* 


■, and thus suboptimal for (P2'). 
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